Jet quenching probe of gluon plasma formation 
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Abstract 

We study, in a pQCD calculation augmented by nuclear effects, the jet en- 
^ ' ergy loss needed to reproduce the vr" spectra in Au+Au collisions at large pT, 

measured by PHENIX at RHIC. The transverse width of the parton momen- 
Qh! tum distributions {intrinsic kx) is used phenomenologically to obtain a reliable 

baseline pp result. Jet quenching is applied to the nuclear spectra (including 
shadowing and multiscattering) to fit the data. 



1 Introduction 



The average energy loss of jets in central Au + Au collisions at RHIC can be used to 
estimate the density of gluons in the medium. To determine the effect of jet energy loss 
from the data, one needs a reliable background calculation of pion and kaon spectra 
subject to jet quenching. Here we provide the required baseline pQCD calculation 
augmented by the nuclear effects of multiscattering and shadowing. 

2 Parton model and pQCD with intrinsic trans- 
verse momentum 

The invariant cross section for the production of hadron h in a. pp collision is described 
in the pQCD-improved parton model on the basis of the factorization theorem as a 

*Presented at the International Europhysics Conference on High Energy Physics (HEP2001), Bu- 
dapest, Hungary, 12-18 July 2001. JHEP, PRHEP-hep2001/242. 
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convolution|jl[]: 



-C/fe-p— = 2^ / dXa,b dZc fa/p{Xa, kxa, Q > fb/p{Xb, krb, Q ) -TJ ' ^ sd{s+t+U) , 

^ P abed •' 

(1) 

where fa,b/p{x, k^, Q"^) are the parton distribution functions (PDFs) for the colhding 
partons a and h in the interacting protons as functions of momentum fraction x, at scale 
Q, incorporating a Gaussian transverse momentum distribution. Having convinced 
ourselves that it makes sense to use a pT-independent width at a given energy |0, 
we treat the width of the transverse momentum distribution, {k'j^)pp, as an energy- 
dependent parameter (also see e.g. [§]). 

In eq. (|T]), dcr/dt is the hard scattering cross section of the partonic subprocess 
ab — >• cd, and Dh/dzc, Q'^), the fragmentation function (FF), gives the probability for 
parton c to fragment into hadron h with momentum fraction Zc at scale Q'. The scales 
are fixed in the present work as Q = pt/2 and Q' = pt/{2zc). We use leading order 
(LO) partonic cross sections, together with LO PDFs (GRV)[|] and FFs (KKP)||. 
This ensures the consistency of the calculation. Approximate NLO treatments via a 
Pt and energy independent "K factor" are also possible [^, |^. However, the dependence 
of the K factor on px and y/s can not be neglected if a wide range of these variables is 
to be addressed[§]. In addition, we prefer to keep the effects of the K factor separate 
from those of the intrinsic transverse momentum. 

Fig. 1 displays the extracted values of the width {k^)pp from hard pion production 
data at energies below -^i < 60 GeV and from = {h~^ + h~)/2 data at higher 
energies. The solid curve is drawn to guide the eye, while the dashed lines delineate 
an estimated band of increasing uncertainty with increasing ^/s. We use this band to 
read off the {k^)pp values to be applied at RHIC energies. The solid bars represent 
estimates at i/i = 130 and 200 GeV, respectively. We use {k'^)pp= 2 GeV^ in our 
pQCD calculations analyzing recent RHIC results at ^/s = 130 AGeV. 



3 Nuclear shadowing and multiscattering 

We approximately include the modification of the PDFs inside nuclei ("shadowing") 
and the isospin asymmetry of heavy nuclei into the nuclear PDFs via the average 
nuclear dependence and a scale independent shadowing function Sa/A{x) adopted from 



fa/A{x,Q'^) = S, 



a/A 



X] 



jfa/p{x,Q')+[l 



-J ) fa/n{x,Q 



(2) 



where fa/n{x, Q"^) is the PDF for the neutron. 

The enhancement of the width of the transverse momentum distribution of par- 
tons due to interactions in the medium ("multiscattering") is treated in the Glauber 
framework ||10||. We hold multiscattering accountable for the Cronin effect in pA colli- 



sions, and use pA data to parameterize the width enhancement. The enhanced width 
is written as 



{krp)pj^ 



{kl)pp + C-iuA{b)-l) , 



(3) 



where {k^)pp is the width of the parton transverse momentum distribution in pp col- 
hsions, z/^(&) — 1 is the number of effective nucleon-nucleon (NN) colhsions at impact 
parameter 6, which impart an average transverse momentum squared C. In pA reac- 
tions, where one of the hard-coUiding partons originates in a nucleon with additional 
NN colhsions, we use the pp width from Fig. 1 for one of the colliding partons and the 
width for the other: 

dp J (xp 

where tA{b) = f p{b, z) is the nuclear thickness function normalized as / <^^htA{h) = 
A. The thickness function determines the collision number z/^(6) as VAip) = aNNiAip), 
where ctattv is the inelastic nucleon-nucleon cross section. We obtain the best fit to the 
Cronin data if we maximize the number of effective NN collisions at a given impact 
parameter at — l)max = 3, and take the value of the parameter C = 0.4 GeV^. 



4 Jet quenching 

In AA reactions, where both hard- colliding partons originate in nucleons with addi- 
tional semi-hard collisions, we use the width (0) for both initial partons. Thus, 



da: 



AB 



d> 



d%d\ tAir) tBi\b-f\)Ef, 



d^p 



(5) 



with all parameters as fixed earlier. Applying this model to WA80|TT| and WA98[0 
central production data at SPS energies, we find a reasonable agreement beyond 
a transverse momentum of pt > 2 - 2.5 GeV, with possibly an overestimation of the 
data by the model in the heaviest system (Pb + Pb) by upto 40% [^. This may be taken 
as a hint that an additional mechanism is at work in the nuclear medium, which acts 
to reduce the calculated cross sections. As we will see next, the discrepancy becomes 
more serious at RHIC energies. Since this effect can also be looked upon as a shift of 
the spectra to lower px, i.e. jet energy loss, we study jet quenching |T3|, as a 
potential explanation. 

Assuming a static gluon plasma, we use here the average energy loss of gluon jets 
to estimate the effect on the cross sections. The fragmentation part of (|l|) is modified 
as 



r2\ 



i2\ 



TTZt. 



TTZt. 



(6) 



where the down-shifted (quenched) momentum p* = Pc — ^E results in a momen- 
tum fraction z* = z^ (1 — IS.E jp^. For the energy shift IS.E we use the "thin plasma" 
(GLV) approximation|T^, appropriate for small values of the opacity n = L/\, where L 
is the average thickness of the static medium and A denotes the gluon mean free path. 
One important feature distinguishing the GLV energy loss from the "thick plasma" 



(BDMS) limit [14] is the approximately constant fractional energy loss (AE/E) of glu- 



ons as a function of jet energy E in the range 3 < < 10 GeV[|T5 . 



In Fig. 2. the results of this calculation are compared to recent 7r° spectra measured 
by PHENIX |ll6| , p!7| for values of the opacity < n < 4. It can be seen that, while 
the semi-peripheral spectrum displays only a weak (but non-vanishing) influence of 
jet-quenching characterized by n ~ 1 - 1.5, in the central case a much larger effective 
opacity of fi ^ 3 - 4 is required for agreement between the data and pQCD calculations 



T8[| . It should be noted that fluctuations in jet energy loss reduce the suppression, but 
their net effect can be absorbed in a rescaling of the opacity n. First results on an 
approximate treatment of fluctuations and on a dynamical calculation are reported in 
the next talk[P^. 

The average opacity n can be related to the density of the plasma using the expres- 
sion of the transport coefficient in terms of dNg/dy [|15|]. Using this correspondence 
implies a very high initial gluon density for fi ^ 3 in a static plasma. Taking into 
account the rapid transverse collective expansion of the plasma leads to somewhat re- 
duced gluon densities, dNg/dy ^ 500, and to results on the collective flow consistent 
with the preliminary data|2^. By any reasonable estimate these gluon densities indi- 
cate that strongly-interacting matter has reached the deconfined region in the analyzed 
nuclear collisions. 



5 Summary and conclusions 

Jet quenching is needed on the background of a pQCD calculation augmented by 
the nuclear effects of shadowing and multiscattering to bring the calculated results in 
agreement with the PHENIX data on 7r° spectra from Au + Au at ^/s = 130 GeV. The 
amount of required jet energy loss indicates a dense gluon plasma, well in the deconfined 
region. Jet quenching is a powerful diagnostic tool to study strongly-interacting matter 
in nuclear collisions, which will play an important role in analyzing the 2001 RHIC data 
at = 200 GeV, and beyond. 
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Fig. 1. The best fit values of {k^)pp in pp nX and pp — > /i^X reactions. Where 
large error bars would overlap at the same energy, one of the points has been shifted 
slightly for better visibility. The band is drawn to guide the eye. 
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Fig. 2. Semi-peripheral and central vr*^ production in Au + Am collisions at ^/s 



130 GeV. Data are from Ref. IjT^- Calculation with jet quenching at opacities 
L/X = 0,1,2,3,4. Dashed lines are pp —>■ vr^ and Glauber-scaled estimates for 
Au-FAu. 



